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Abstract 

For two gaussian states with given correlation matrices, in order that relative entropy between them is 
] practically calculable, I in this paper describe the ways of transforming the correlation matrix to matrix in the 

exponential density operator. Gaussian relative entropy of entanglement is proposed as the minimal relative 
entropy of the gaussian state with respect to separable gaussian state set. I prove that gaussian relative 
entropy of entanglement achieves when the separable gaussian state is at the border of separable gaussian 
state set and inseparable gaussian state set. For two mode gaussian states, the calculation of gaussian 
relative entropy of entanglement is greatly simplified from searching for a matrix with 10 undetermined 
parameters to 3 variables. The two mode gaussian states are classified as four types, numerical evidence 
strongly suggests that gaussian relative entropy of entanglement for each type is realized by the separable 
state within the same type. For symmetric gaussian state it is strictly proved that it is achieved by symmetric 
gaussian state. 

PACS number:03.67.Mn, 03.65. Ud 
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1 Introduction 

CN 

Quantum relative entropy function has many applications in the problems of classical and quantum information 
transfer and quantum data compression P^. The relative entropy has a natural interpretation in terms of the 
statistical distinguishability of quantum states; closely related to this is the picture of relative entropy as a 



o 

[ distance measure between density operators. Based on the relative entropy, a nature measure of entanglement 



X 



called the relative entropy of entanglement was proposed. This entanglement measure is intimately related to 
the entanglement of distillation by providing an upper bound for it. It tells us that the amount of entanglement 
in the state with its distance from the disentangled set of states. In statistical terms, the more entangled a state 
is the more it is distinguishable from a disentangled stated- However, except for some special situations [3], 
such an entanglement measure is usually very difficult to be calculated for mixed state. For continuous variable 
system, it is shown that the relative entropy of entanglement is actually trace-norm continuous and hence well- 
defined even in this infinite-dimensional context 0]. Due to the fact that gaussian state has the merit that the 
logarithmic of the state operator is in the quadrature form of canonical operators, the distance of the gaussian 
state to gaussian state measured by the relative entropy (gaussian relative entropy or GRE) was considered [5] . 
In which the second gaussian state was specified by its exponential operator matrix (EM), not by the usual 
correlation matrix (CM). Until now, an explicit transform from the CM to EM has not been available, make 
the calculation of the relative entropy between gaussian states in fact impossible. Moreover, when using the 
gaussian relative entropy as entanglement measure, the second state should be separable. Now all the separable 
criterions are not in the form of EM. So the problems of separable criterion in the form of EM or the explicit 
transform of CM to EM need to be addressed. In this paper, I will give the explicit transform of CM to EM, 
propose the gaussian relative entropy of entanglement (GREE) and give the method of how to calculate it. The 
paper is organized as follow: In section 2, the transform of CM to EM for q — p decorrelation gaussian state is 
given with the help of symplectic transformation, further more, the direct transforms of any CM to EM and EM 
to CM are also solved due to the commutation relation of matrices. Section 3 deals with GREE with emphasis 
on the proof of theorem 2 which states that GREE will achieve by gaussian state at border of separable and 
inseparable gaussian state sets. In section 4, I concentrate on the simplification of GREE for two mode (1 x 1) 
gaussian state system. In section 5, the lxl gaussian states are classified as four types, their GREEs are 
discussed. In section 6, conclusions and discussions are addressed. 



1 



2 Matrix in the exponential density operator 



To characterize a gaussian state, we have several equivalent means, among them are: quantum characteristic 
function specified by first (irrelative to the entanglement problem) and second moments which are also called 
means and CM, density operator in exponential form specified by a matrix M (exponential operator matrix or 
EM), density operator in exponential form of ordered operators specified by another matrix (ordered exponential 
operator matrix or OEM). The separability of a gaussian state was obtained with CM 00 0, also with OEM 
[Hj. The transform of CM to OEM and vice versa are quite directly by the integral within ordered product 
of operators. The transform of EM to OEM is also available but involved with a calculation of exponential 
of matrix^U] . Scheel [5] derived a relation of EM and CM of gaussian state with generation and annihilation 
operators. Following the way I now derive their relation with canonical operators. 
Gaussian quantum state can be given by the density operator 

exp[-±F T MF] 

(1) 



Tr{exp[-±F T MF}}' 



where M is a real symmetric matrix, and F = [q\, ■ ■ ■ q n ;p\, ■ ■ ■ p n ] , with qj,Pj are the canonical operators. 
In order to relate the matrix M to the correlation matrix (CM) a, a unitary transformation 

F = UFU- 1 = SF (2) 

is needed. The matrix S produces a symplectic transformation on the canonical operators F. To preserve the 
commutation relations of the canonical operators, S should satisfy the symplectic condition SAS T — A. Where 



/„ 

-In 



(3) 



The matrix S is chosen such that it diagonalizes M , hence S T MS = M (with M being diagonal). Then the 
characteristic function of the density operator (JIJ is 



Tr[pexp(iF T z)] = Tr[U pU~ l U exp(zF T ' z)U] (4) 
= eM~z T S(±coth±M)S T z], 

and the CM is a = S(h coth |M)S ,T = SaS T . In the derivation the expression for characteristic function of a 
thermal state has been used. 

For a given CM, how to find such a sympleptic transformation is the topics of this section. To find the S 
matrix, let us consider A _1 a instead of a, then 

5 T A-V(5 T )- 1 = A- 1 5. (5) 

The eigenvalues of A _1 a come in pairs ±«7j, the matrix uncertainty relation requires jj > | |13| . where jj 
are called symplectic eigenvalues of a then a = diag{ji, • • • , j n , 71, • • • , 7„}, and M — diag{M\, ■ ■ ■ , M n , 
Mi, • • • , M n }, with Mj — log 2^~j are symplectic eigenvalues of M. Let 4" be the eigenvector of A _1 5 and ^ 
be the eigenvector of A _1 a, then 

* = (5 T )" 1 $. (6) 

So that S matrix can obtained from the eigenvectors. The S matrix so obtained is not totally determined, yet 
the symplectic condition SAS T = A should be verified. 

If the CM is in its q — p decorrelation form, that is a — a q © a p , then the details of S can be further worked 
out. The q—p decorrelation CMs are quite general in concerning with entanglement for lxl gaussian states. For 
these states, the CMs can always be transformed into the form of a q © a p by local operations [S| [7j. For multi- 
mode bipartite gaussian state, It is not known if the CM can be transformed to the form of a q © a p or not just 
by local operations. Let us start with the CM of the form a q ©a p , that is to say the correlation between position 
and momentum of each mode and inter-modes have already dissolved, so that one just considers a symplectic 
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transformation of form S — S q ® S p . The relation SAS T = A then will requires that S q S p — S p S q — I. And S 
will take the form of S = S q ® (S'J)" 1 , and (S^" 1 = (S'J)" 1 ©S , ? . S^a^)- 1 = [S^a^^)" 1 ] ® [Sja p S g ] - 
5. Let ty T = (\E r ^,'&T) an( j ^j,t _ (i&T ^,T^ then the eigenequation A~ 1 a<I' = ±17^ (where 7 is one of the 
symplectic eigenvalues of a) will be —a p ^ p = ±ijfy q and a q ^ q — ±17^. And one gets 

a p a q ^ q = 7 2 * 9 , a 9 a p * p = 7 2 4y (7) 

For the eigenvalue ijj, suppose the eigenvector of the q part be = [cji t • • • , Cj n ]. Since a p a q ^ q = 7 2 * 9 , the 
eigenvector ^> q can always be chosen to be real. Then = — z[cji, • • • , Cj n ]a g /7j-. The similar equation for ^> 
will at last give the result of 

$J = —= [0, • • • , 0, 1, 0, • • • , 0; 0, • • • , 0, -i, 0, • • • , 0] for eigenvalue ij J: (8) 
v2 

where the nonzero elements in the eigenvectors are at the positions of j and n + j. By Eq.JSJl one gets 
* T (S ,T )- 1 $ = and $+(S' T )- 1 $ = then 



q >kj 



(S, 
{S q )kj 



Cjk 

£■ 



One has to verify the consistence of Eq.® and Ea. (|lt)(l . that is S q l S q = I n . The vector 
factor and its length left to be determined. To determine the phase factor of \P 



qj , one simply chooses Cjj 



(9) 
(10) 

has yet a phase 
to be 



positive. The length of the vector ^ q j can be chosen such that 



kl 



CjkCji{a q )i k /jj 



(11) 



' mk 



(S q ) kj = for m ^ j, one has 



.quip 



-a q a p )^ qm = 



pj ((a p a q ) - a q a p )$fq m = (7? - 7^) = 0. Hence, for 7 m ^ 7/., one gets *+-f gm = 0, and the proof 



So Ek {S- 1 )^ (S q ) kj = 1. While in proving £ fe (S" 1 

of the consistence of Eq. © and Eq. (|10|l is completed 

Hence the symplectic transformation can be constructed directly from the eigenvectors of a p a q or equiv- 
alently a q a p . As an example, let us consider the lxl symmetric gaussian state. The q and p parts of CM 
are 



1 

a "= 2 

The symplectic transformation will be 



rn 



rn 



where s\ = 



S2 = 



S 



7n — k q 
m+k % 



1 


Si s 2 


1 


' 1/si l/s 2 


71 


si -s 2 




. 1/Sl -1/82 _ 



And one obtains the M-matrix for the state 



M = {S T )- 1 \og 



25+1 
25-1' 



S~ 



(12) 



(13) 



(14) 



where 5 = a q © 5 P , 5 9 = 5 P = ^diag{y' \m + k q ) (m — fc p ), ^/ (m — fc g ) (m + 

The above transform is limited to the q — p decorrelation matrices. For a more general CM to EM transform, 
although the symplectic transformation is also available up to the length and phase factor of each eigenvector of 
A _1 a, the symplectic condition SAS T = A is not easily verified. So let us turn to a direct way of transforming 
CM to EM as well as EM to CM. This is due to the fact that a = SaS T and S T MS = M. The late can be 
rewritten as M -1 = S M~ 1 S T . So that a and M -1 can be simultaneously symplectic diagonalized. A _1 a and 
A~ 1 M~ 1 will have the common eigenfunctions. And they commutate with each other, [A _1 a, A -1 M _1 ] = 0. 
Hence 



MaA~ 



A aM = 0. 



(15) 



Together with a = \ coth hM, one can transform a to M or M to a. 
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3 GREE and border state 



Now given a CM, one can transform it into EM. This enables the calculation of relative entropy between two 
gaussian states to be practically possible. The relative entropy of a gaussian state p with respect to another 
gaussian state a is defined as 

S(p\\a) =Trp(logp- log a). (16) 
The normalization factor of the state cr is 

c-n^St-n-j-j. (i7) 

3 = 1 3 = 1 - 4 

Hence 



- Trp log cr = - log c + —TrpF T M a F 

= -\ogc+^Tr[a p - l -^jM a 

= -\ogc+^Tra p M a , (18) 

where the operator trace TrFpF T = [a p - |A) [TT] and the fact that TrAM a = TrA(S T y 1 M a S- 1 = 

TrS~ 1 A{S T )~ 1 M a = TrAM a = have been used. The trace in the first equality is two fold, operator trace 
and matrix trace. 

The relative entropy of entanglement was defined as the minimization of the relative entropy of a state with 
respect to all separable state: Er{p) = min^c S(p \\a), where D is the set of separable state. If its subset Dq 
of all gaussian state is used instead of the set D itself , then the GREE for a state can be defined as : 



E GR (p) = min S(p\\a). (19) 
o-eDa 

I will prove that the separable set can be further restricted to the border separable set. For completeness I will 
start theorem 1. 

Theorem 1: The relative entropy of entanglement is obtained when the separable state is at the border of the 
set of separable states and the set of inseparable states. 

Proof: The relative entropy is jointly convex in its arguments |12|. That is, if pi, p2, o~\ and o~2 are density 
operators, and p\ and pi are non-negative numbers that sum to unity (i.e., probabilities), then S(p\\a) < 
PiS (pi ||<7i ) +P2S (p2 1 1 o"2 ) , where p = p\p\ + P2P2, and a = p\U\ + P2°~2- Joint convexity automatically implies 
convexity in each argument, so that 

S{p\\a)< Pl S (p Ha! ) + P2 S (p \\0-2). (20) 

and for < x < 1, one has S (p || (1 — x) p + xa ) < (1 — x)S (p \\p) + xS (p \\a) — xS (p \\a) < S (p ||er) . Hence 
for a separable state a that is not at the border one can find a new separable state with less relative entropy 
until the new separable state is at the border. 

Theorem 2: The gaussian relative entropy of entanglement for gaussian state is obtained when the gaussian 
separable state is at the border of the set of separable states and the set of inseparable states. 

Proof: The idea is like this: for any given separable gaussian state 00, one needs to find a line to connect ao 
and the inseparable gaussian state p with every point in the line is a gaussian state which is denoted as cr. In 
the line, between the separable state ao and inseparable state p, there should be a border gaussian state. I will 
find such a line by continuously change the state a in the fashion that the relative entropy of p with respect 
to a decreases monotonically. If the process of decreasing of relative entropy does not stop, then the relative 
entropy will go to its minimum value. Because S (p \\a) > and with equality iff p — cr, so a will eventually 
reach p. In the following I will mainly prove that the decreasing process would not stop if cr ^ p. 

Now unitary operations leave S (p\\a) invariant, i.e. S(p\\a) = S (UpU + \\UaU + ) . This reflects the fact 
that 

Tra p M a = Tra p (Sj)" 1 M a S~ l = TrS^ctp (Sj)' 1 M a . (21) 
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Denote (3 = S a 1 a p (S'J) . The relative entropy will be 

n ( M \ 1 ™ ~ 
S(p\\a) = Trplogp- £ log 2 sinh -f- + - £ + M„, (22) 

Because /3 is also a CM of some gaussian state, the uncertainty relation requires (3— |A > 0. Hence Pjj(3 n +j,n+j > 
\, and \ ((3 j j + f3 n +j,n+j) > y/ /3jjf3 n+j , n +j > \- Denote (3 oj = \ (j3jj + (3 n+jtn+j ). The partial derivatives of 
the relative entropy are 

dS(p\\a) 1 ln n . 1 , M CT j 

= Pjj-lrt- (23) 

Because M CT j is a monotonically decreasing function of j a j , the partial derivative of the relative entropy with 
respect to 7 CT j is positive iff 7^ > The line designed to connect Co and p is like so: first let us fix S a to 
S ao , then fix all j a j to 7 CTo j except j a i, if 70-1 > then decrease 7 CT i until it is equal to the relative 
entropy decreases monotonically. if 7 CT i < /3 n , then increase 7 CT i until it is equal to /3 11} the relative entropy 
also decreases monotonically. Now the state is with 70 1 = /3 n and all other 70-/ = j ao j- Then let us make 7^2 
monotonically vary to /3 2 2 while keeping the other j a j fixed. At last all j a j — (3^, and in every step the relative 
entropy decreases monotonically. Now the state o is with its all symplectic eigenvalues j a j = /3 • • but with So- 
still being fixed to S ao . The relative entropy then will be 

S(p\\a) = Trp logp + £ 9 ±) , (24) 

where g (x) = (x + 1) log (a; + 1) — xlogx is the bosonic entropy function, which is a monotonically increase 
function of its argument, but its derivative d9 }^ decreases with x increases. 

The next step of stretching the line is to change S a gradually in order that S (p \\a )decreases further. 
Now (3 = S^ 1 SpdtpSj (Sj) , let us apply infinitive small symplectic transform to (3 to change the state a 
continuously. The infinitive small symplectic transforms will accumulate some finite symplectic transforms, 
which are the following six kinds: (i) Local rotations which keep (3jj invariant; (ii)Local squeezings which can 
be used to decrease /3 - • to y/f3jj/3j +n ,j+n ; (hi) The first kind two mode rotations. For modes i and j, if the 
two mode CM (3 su b (submatrix of (3) is arranged according to the order of canonical operators [qi, qj]pi 7 pj], the 
rotation will be 9 (0) © 6 (9), where 



e(0) 



cos 9 sin 9 
— sin 9 cos 9 



(25) 



Before the rotation is applied, (3 su b has already been prepared (local squeezed) to the form of equal diagonal 
elements within each mode, that is (3 U = (3u — f3i+ n ,i+n and the same for mode j. The inter- mode rotation keeps 
(3 U + f3jj that is the trace of the two mode CM (3 su b invariant. The only way to decrease the relative entropy 
is to enlarge the difference between (3 U and f3jj. This is because that the bigger one say (3 U increases some 
amount, the smaller one f3^ will decrease the same amount, but total relative entropy will decrease due to the 
monotonically decreasing property of the derivative of bosonic entropy function dg j x>> . The distance between (3 U 



dx 

and fljj can be enlarged at most to 



3) 



P u - Pjjl 2 + {fiij + (3i+ n ,.j+n) 2 by proper rotation. After 

such rotation, we have (3^ — —0 i+n j +n , the off diagonal elements of q part and p part are asymmetrized ;(iv) 
The first kind two mode squeezings R(r) © R (— r), where 



R(r) 



cosh r sinh r 
sinhr coshr 



(26) 



After successive applying of (ii) and (iii) for rounds, (numeric results indicate the a few rounds will do), (3 su b 
will have the form with (3u — (3i+ nt i+ n , Pjj — f3j+ n .j+n,l3ij = —(3i+ n ..j+n, The two mode squeezing then will be 



■5 



used to diagonalize this part of (3 su b matrix. The squeezing decreases f3 u and /?. ,,- with the same amount which 



■\f\Pii ~ Jjjf + (ftj'+n - Pi+n,.j) 2 - And 



33 

can be at most + /3jj) — yj \{f3 u + Pjj) 2 — Pfj ■ ( v ) The second kind two mode rotations which rotate 

qi,Pj pair and simultaneously Pi,qj pair. The rotation is similar to the first kind two mode rotation but with 

distance between j3 u and (3^ can be enlarged at most to (3 U — (3j 

one has j3 i - +Jl = fi i+n j after the such rotation, (vi) The second kind two mode squeezings which squeeze qi,Pj 
pair and simultaneously pi, qj pair. After successive applying of (ii) and (v) for rounds, f3 su b will have the form 
with f3u = fii+r^i+m f3jj = flj+nj+m f3i.j+ n = Pi+n,.j, The two mode squeezing then will be used to diagonalize 
this part of f3 su b matrix. The squeezing decreases (3 H and (3^ with the same amount which can be at most 

\{Pu + Pjj) ~ \J jWu + Pjj) 2 ~ 01 ' j+n- These six symplectic transforms are classified as three group: (i); (ii)- 
(iii)-(iv); (ii)-(v)-(vi). Each group aims at diagonalized 4 of the off-diagonal elements of the two mode CM while 
decreasing the relative entropy except the group (i) which keeps the relative entropy. By successively apply the 
three groups the relative entropy will decrease step by step before it is diagonalized. Then the whole procedure 
of diagonalizing is applied to all other pairs of modes round and round. Before (3 is totally diagonalized, a way 
can always be found to decrease the relative entropy. The totally diagonalized (3 (denoted as (3) is exactly 5. 
This is because the original (3 = S~ 1 a p (S^) has the same symplectic eigenvalues with a p , so (3 and a may 
only differ by the interchange of mode, However the relative entropy is according to Ea. H24|) . This is only 
possible when a — p. So that a at last reaches p. Maybe the gradually changing a passes the border of separable 
and inseparable sets several times, but this does not matter, the last passing meets the requirement of GREE. 
And theorem 2 is proved. 

In finding the minimization in the border state set, one has another question, that is, if displacement decrease 
the relative entropy or not? The answer is negative. The operator is exp (iF T z) , where z is a real 2n vector. 
Since 

exp (iF T z) crexp (-iz T F) = cexp[-i (F + Az f M a (F + Az)} (27) 



1 T 
Trp log a = - logc + -Trp (F + Az) 1 M a (F + Az) 

- - logc + ^Tra p M a + ^Tr (Azf M a Az 



(28) 



While M a is positive definite (I will elucidate it at section 5), so the last term is not less than 0. The displacement 
can not decrease the relative entropy. 

To find the GREE of state p is to find the M a matrix of a border state a such that the relative entropy 
reaches its minimum. 



4 GREE of 1 x 1 gaussian state system 

Now let us turn to the lxl gaussian state system. The general case of relative entropy is that a p is in its 
standard form but M a is not. It is no need to require that they are all in the most general form, because by the 
unitary invariant of relative entropy, at least one of the matrices a p and M a can be converted to any possible 
form. Since any possible M a matrix can be simplified to its standard form by local operations^, so such a 
M CT can be generated from the standard form M as . For lxl system, suppose a p takes the standard form, 





Oti 


a 2 




a\ 


Q!4 




a.2 


a 3 




Qf4 


a 3 



(29) 



M as takes the same form but with elements M S i (i = 1, • • • ,4) respectively. The local operations are firstly 
a local rotation L\ with angles and Obi for the two modes respectively, then a local squeezing L 2 = 
diag{exp (ta) , exp (tb) , exp (—ta) , exp (—tb)}, then another local rotation L3 with angles 9 ai and 9bi for the 
two modes respectively. The standard form of M as is modified to M a = L\lfL\M as L\L2L?,. Local operations 
will leave the normalization factor c unchanged, so one just needs to consider Tra p M a term, one gets 



G 



Tra p M a = 2aiM sl cosh 2t a + 2a 3 M s3 cosh 2t b 

+ cos (9 A i + B i - A2 - B2 ) (a 2 - a 4 ) (M s2 - M sA ) sinhr^ sinhT B 
+ cos (0 A \ + Obi + A2 + B2 ) {ot 2 - a 4 ) {M s2 - M s4 ) coshr^ coshr B 
+ cos (0 A1 + B i - A2 + B2 ) (a 2 - a 4 ) (M s2 - M s4 ) sinhr^ sinhT B 
+ cos (0 A1 + B i + A2 - B2 ) (a 2 + a 4 ) (M s2 - M s4 ) coshr^ sinhr B 
+ cos (0 A1 - B i - A2 - B2 ) (a 2 - a 4 ) (M s2 + M s4 ) sinhr^ coshr B 
+ cos {6 ax - B \ + A2 + B2 ) (a 2 - a 4 ) (M s2 + M s4 ) coshr^ sinhr B 
+ cos (9 A x - B i - A2 + B2 ) (a 2 + a 4 ) (M s2 + M s4 ) sinhr^ sinhrs 
+ cos {6 ax - Obi + A2 - B2 ) (a 2 + a 4 ) {M s2 + M s4 ) cosh r A cosh t b . (30) 

Clearly, in order that Tra p M a is minimized, the local rotations should be arranged in such a way that all the 
cos-factors are ±1, then 

Tra p M a — 2a\M s4 cosh 2t a + 2azM s3 cosh 2t b 

- [\{a 2 - a 4 ) (M s2 - M s4 )\ + \{a 2 + a 4 ) (M s2 + M s4 )\] cosh {\t a \ + \t b \) 

- [\(a 2 - a 4 ) (M s2 + M s4 )\ + \(a 2 + a 4 ) {M s2 - M s4 )\] sinh (\t A \ + \t b \) ■ (31) 

Without lose of generality, let a 2 > 0,a 4 < and a 2 > —a 4l and the M a matrix now has the form of 
L' 2 M' as L' 2 with L' 2 = diag{exp\r A \ ,exp\r B \ ,exp (— \t a \) ,exp (— \t b \)}. M' as only differs from M as by off 
diagonal elements. For simplification of the notations, denote M' asi as Mi, then M 1( - 3 ) = M sl ( 3 ) and M 2 ( 4 ) = 
-\(\M s2 + M s4 \ ±\M s2 -M s4 \). 

The problem now is to determine the elements of the matrix M a of a border state a. The local squeezing 
L' 2 can be rewritten as the product of two local squeezings Y and X, Now M a = Y{y)X{x)M' as X{x)Y{y) 1 with 
Y(y) — diag{y/y, y/y, y/y^ 1 , y/y 1 } and X(x) = diag{y/x, \fx , yfx , yfx\ are symplectical transformations. 
After minimization of Tra p M a with respect to y, one has 

-TratpMv = yJ(axMxx + a^hx' 1 + 2a 2 M 2 )(a 1 M 1 x- 1 + a 3 M 3 x + 2a 4 M 4 ). (32) 

The further minimization of \TrotpMa with respect to x will lead to an algebra equation of x up to power 
of 4. Although after this round of minimization, ^TrQ. p M a can not easily be expressed, but in principle it is 
possible to analytically express it as a function of the four Mi. The GREE problem of gaussian state p now is 
the minimization of — logc+ ^Tra p M a where the whole function now has only the four Mi as its variables. 
The four Mi would fulfill the border state condition, so there are only 3 of them left to be determined in further 
minimization of the relative entropy. Consider the border state condition for state characterized by its standard 
form M' as , that is by the four Mi, suppose the corresponding CM is a = a q © a p , the condition for the state to 
be a border state is [7]: 

4dct (a q a p ) = Tr (a q a p ) + 2 (\cxc 2 \ - cic 2 ) - -. (33) 
where ci, c 2 are off diagonal elements of a q and a p respectively. By Eqs. 0, one has det (a q a p ) — 7^47^ and 

Tr(a q a p ) =7a + 7bi ( 34 ) 

where 7,4, 7b are the symplectic eigenvalue of a for the two modes respectively, and jj — 5 coth (j = A, B). 

While Mj are the symplectic eigenvalues of M' as , what left is to determine c\c 2 . The commutation relation 
[A _1 a, A _1 A/ _1 ] = now takes the form M p a p = a q M q or equivalently a p M p — M q a q . This enables all the 
other elements of a to be expressed as linear combination of ci and c 2 . Since det (a q a p ) — det (a g M~ 1 a 9 M 9 ) = 
[det (a q )] 2 det (M q ) / det (M p ) , so that 



det (a q ) = lAlB JAet (M p ) / det (M q ). (35) 
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Either Tr (a q a p ) or det (a q ) is a linear combination of cf , c| and C1C2. By combining Ea. lfMjl and Eq.(J23), one 
arrives at a quadratic equation about cic 2 . So that C1C2 can be expressed with Mj so does the border state 
condition of Eq. . 

In this section, I reduced the EM of the destination state from 10 parameters to 3 (4 parameters and 1 
restriction, strictly speaking). In the next section I will elucidate that it is really 3. These 3 parameters are 
left for numeric calculation, because minimization function at this step is too complicated to be dealt with 
analytically. 

5 Classification of 1 x 1 gaussian states 

In this section, I will aim at constructing the six parameter M a , the footnote a will be omitted when it is not 
confusing. The CM a should satisfy uncertainty relation a — > 0. The standard form of the correlation 
matrix of lxl gaussian state have four parameters. Uncertainty relation adds some restrictions among the four 
parameters. So that the parameters are not freely chosen, otherwise the state may not be physical. Needlessly 
to say, EM M is less restricted than CM a. If all symplectic eigenvalues of M isrpositive, the state should be 
physical, because the density operator will be in the form of a ~ t/cxp(— \ J2j ^jilj + Pj))U + = AA + , with 
A = U exp(— Mj(g 2 + p 2 )), so that a is positive definite. For standard form M of 1 x 1 gaussian state, 
it is easy to check that M should be positive definite. But as elucidated in the former section, the separable 
criterion is quite complicate expressed with M. In this section, I will seek free parameter representations which 
are simple both in the uncertainty relation and separable criterion. 
Given a standard form of a = a q © a p , with 





a 


Cl 




a — c 2 


a q = 




b 




— c 2 b 



(36) 



where without lose of generality c\ > and — c 2 < are supposed ( when c±(— c 2 ) > 0, the state is def- 
initely separable, so that will not be the border state of interesting). Now let us seek an operational way 
to symplectically diagonalize it. This is accomplished by first applying local squeezing X, then a two mode 
squeezing or a two mode rotation according to different structure of a. The CM will be transformed to (i) 
R(r)®R (-r) X{x)aX(x)R (r) © R (-r) or (ii) 9 (0) © 6 (g) X(x) aX(x)Q (-0) © 6 (-0) . In case (i), x is so 

chosen such that (ax + bx^ 1 ) / c\ = (ax^ 1 + bx)/c2, x — ^/ °^~^ , then a two mode squeezing with tanh(2r) 

= —2y/{ac\ — 6c 2 )(ac 2 — bc\)/ |a 2 — 6 2 | will diagonalize the properly local squeezed CM. Clearly the existence 
of x requires that f + f > + f^- In ca se (ii), x is so chosen such that {ax — bx^ 1 )/^ = -(as -1 — 6x)/c 2 , 

x = ij fee i- act ' t nen a two mode rotation with tan(20) = 2^J (ac\ — 6c 2 )(&ci — ac 2 )/(a 2 — b 2 ) ■ sign(bci — ac 2 ) 

will diagonalized the CM. The existence of x requires that f + f < + f^- After the diagonalization of these 
two cases, further local squeezing of X(x') and Y(y) will be applied to transform their diagonal elements into 
symplectic eigenvalues. Certainly there are the third case of f + £ = f^ + f^ and the fourth case of a = b 
which is the case of symmetric gaussian states. So, according to | + | being more than or less than or equal to 
f~ + ^r, the state is classified as type (i), type (ii), type (iii). And if a — b, the state will be type (iv) state. The 
quantity (| + ^)/(fj + f 2 -) is critical in symplectic diagonalization of the CM. It is a kind of ratio of diagonal 
elements to off diagonal elements of the CM. 

Now let us construct the CM of all four classes. The process is just the reverse of the diagonalization. Let 
us begin with a = diag{jA,jB',^A,jB}, Then X(x'^ 1 ) is applied. Y(y~ 1 ) can be put at the last stage because 
it commutates with all other kind operations of q — p decorrelation type such as the two mode rotation and the 
two mode squeezing and local squeezing X. The next step is to apply R(—r)®R (r) or G (—0) © G (—0) for type 
(i) and type (ii) CMs or states respectively. The separable property is totally determined after this step and the 
successively applications of local squeezing operations X(x~ 1 ) and Y(y^ 1 ) will not affect the separability. So 
when the separability is concerned, the last two local squeezing operations can be omitted. The CM generated 
will be (i):R (-r)ffii? (r) X(x'- 1 )5X(a;'' 1 )i? (-r)ffii? (r) and (ii) G (-0)©G (-0) X{x'- l )aX{x'- l )Q (0)©G (0) . 
The separable criterion for the two cases then will be 

(27a " -\)> sinh 2 (2r)[(x' 2 + x>~ 2 ) lAlB + ( 7 2 + 7 |)] (37) 
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for type (i) and 



(27l " ^)(2 7 I - > sin 2 (20) [(x' 2 + a;- 2 ) 7 a7b - ( 7 a + 7 2 )] 



(38) 



with x' > maxf./^, \/~} or ^ < m i n {\/— ; \/~} f° r type (ii), where the equalities in these two equation 

are for the border states, For border state, one of the parameters, say x', can be easily expressed by the other 
three parameters. The corresponding EMs will be (i):Mj = R(r) © R(-r) X(x')MX(x')R{r) © R(—r) and 
(ii) M n = © (-#) © © (-0) X(x')MX(x')Q (8) © 9 (0). The EMs so generated are usually not in the standard 
form, but this does not matter, after local squeezings X(x) and Y(y) being applied, the six parameter form EMs 
which are the most general of q — p decorrelation M a will be generated. After minimization of the Tra p M a with 
respect to y and x as described in the former section and by using of the border state condition, the relative 
entropy will become a function of 3 variables which are JAi7b,t for type (i) states and 7 a,7b,0 for type (ii) 
states. 

The type (iii) CM with f + | = f~ + f- and type (iv) CM with a = b can be converted to their EMs by 
solving the symplectic transformation matrix S directly as in section 1. The EM of type (iv) state is already 



given in section 1, with the border state condition of (m 
form with only two variables 7a, 7s through 



k aq )(m a — k ap ) = 1, the border EM will be in a 




4-, 



47b ( 4 7a + 1) 



The type (iii) CMs are of two kinds: | = and | = The S q matrices will be 



(39) 



71 



and 



0, 



for the two kinds after applying the border state condition for each, where 5 = ( 7 ^ — j)('y% — \)- 

The numerical results of all four type states indicate that GREE of the state will be achieved by the border 
state of the same type. I calculate the minimization of GRE of a given state with respect to all four type 
border states. Some of the results of type (i) and type (ii) states are displayed in Fig. (1) and Fig. (2). For 
type (iii) and type (iv) states, the border state which achieves the GREE is also type (iii) state or type (iv) 

state respectively. Moreover, it is worth noting that the state realized GREE has the value (f + f ) / f f£ + 

which is very close to but not more than that of the original state p. The numerical results are displayed in 
Fig. (3). This adds evidence to the former numerical conclusion on types. 

One the other hand, for symmetric Gaussian state (the type (iii) state), Based on the numerical observation, 
the conclusion can be drawn to be: for each type of state, GREE will be achieved by the same type of state. 

For type (iii) state ps (symmetric gaussian state) characterized by CM of Eq. (|12|l and further by to, k q ,k p , 
the conclusion that GREE is achieved by symmetric Gaussian state can be strictly proven. Here I will give the 
main idea of proving. The detail will appear elsewhere. Let us start with the 6 parameters^Mo- as was__done_at 
the beginning of this section. This 6 parameter EM can be produced as well as reduced to M a = diag{MA, Mb, 
Ma, Mb} by symplectic transformation S, S T M <7 S — M a , Now the S matrix has 4 parameters which are 
independent of the symplectic eigenvalues. The S matrix takes the form of S — S q © (S^)^ 1 It is known that 
symplectic transformation can always be dissolved to a rotation and squeezing then a successive rotation, that 
is S = R2DR1. with i?i and R2 are rotations and D = diag{d,A, Ab, 1/^a, is the squeezing operation. 

Now for the q — p decorrelation form of S, S = S q © (S'J) -1 , R\ and R2 will be in their simple form of 
i?i = 9 (6) © 9 (8) , and R 2 = 9 (0) © 9 ((p) , By minimizing the relative entropy with respect to 9 and cj) under 
the restriction of a being a border state, one gets 



d fl 

88 { 2 Tra P Ma + XciC2 



d fl 



0, 



(40) 



Where Xjs the^multiplier. In general, the solution is quite complicate and involved with the other four parameters 
d,A,d,B, Ma, Mb- But for the symmetric Gaussian state ps, the solutions are quite simple. They are (i) s'm8 = 0, 
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3.5 




0.2 0.4 0.6 0.8 1 

(2y A -1)/(2 V 1) 



Figure 1: For type (i) gaussian state with sinh(2r) = sinh{2rborder) + 5, a; = 1.1, (2jb — 1)/(27b + 1) = 0.5. 
Solid line is for the searching result of type(i) border states which achieve GREE, dash for type(ii) state reaching 
the minimum of GRE within the type, dotted for type (iii),and dotted dash for type (iv). 
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Figure 2: For type (ii) gaussian state with sinh(26) — 0.5, a; = Xb or d er + 1-5, (2"/b — 1)/(27b + 1) = 0.5. Solid 
line is for searching result of type(ii) border states which achieve GREE, dash for type(i) state reaching the 
minimum of GRE within the type, dotted for type (hi), and dotted dash for type (iv). 
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Figure 3: Solid lines are for original gaussian states whose sinh(2r) = sinh(2rt, or der) + 5, X = 1.1, dotted for the 
border states which achieve GREE, In the left from up to down are for (27s — 1)/(2jb + 1) = 0.2, 0.4, 0.6, 0.8 
respectively 



tan</> = ±l;(ii) cos (9 = 0, tanc/> = ±1. And the other four parameters d,A,d,B, Ma, Mb are not involved. The 
state a then will be the symmetric Gaussian state. 
The GREE will be 

Egr(ps) = Trps log ps + _min {- ^ log(2sinh^-) (41) 



M A ,M B 



3=A,B 



+ + k q )(m - k v )M\ + (m - k q )(m + k p )Ml 

— — Ma Mb 
+ (m — k q )(m — k p )MAMs coth — — coth 

— — Mb 1 
+(m + fc g )(m + k p )MAM B tanh tanh —^-} 2 , 

where Mj = log f-jr-rj , and 7j are symplectic eigenvalues of border state. Here it is easy to prove that no further 
X{x) operation is needed when the border state is prepared in its EM with the form of Eq. (|14(l (but with 
different parameters), that is to say x = 1. 

If k p = k q , the state will be two mode squeezed thermal state^Hl Pst- Then Eq. (|41|) is symmetric for Ma, 
Mb- Clearly the minimum will be achieved at Ma = Mb- So that 

EgrXpst) = TrpsT^ogpsT + inin{-21og(2sinh — -) (42) 

+^Mj4[(m — k q ) coth + (m + k q ) tanh — ^— ]• 



6 Conclusions and Discussions 

Gaussian relative entropy of entanglement is an entanglement measure in its own right. The relative entropy 
between two gaussian states was expressed as correlation matrix of the first state and matrix in the exponential 
density operator of the second state. The mutual transform of q — p fashion CM and EM was derived with 
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symplectic transformation for q — p decorrelation state. A most general transform of CM to EM and vice versa 
was given through commutation relation of the matrices and relation between the symplectic eigenvalues of 
the matrices. I proved that gaussian relative entropy of entanglement achieves when the separable gaussian 
state is at the border of separable and inseparable sets. The displacement or first moments of the second state 
can be ruled out as far as GREE is concerned. For GREE of 1 x 1 gaussian state, the ten parameters EM 
of separable state which minimizes the relative entropy was reduced to three variables EM. Where the matrix 
was decomposed as local operations applied to a standard form of EM. The three variables in EM were left 
for numerical calculation of the minimization. To construct an EM more suitable for the calculation of GREE, 
I classified the standard form CM of 1 x 1 gaussian state into four types according to some kind of ratio of 
diagonal to off diagonal for the first three types and symmetry for the fourth. The numeric evidence on the 
minimization of EMs strongly suggests that GREE for each type of gaussian state will be realized by the state 
within the same type. 

I strictly proved that GREE for symmetric Gaussian state is achieved by symmetric gaussian state. It was 
given as the minimization of a function on the two symplectic eigenvalues of EM. Although the minimization 
equations are easily obtained, but they can not be solved analytically. Further more, for a special kind of the 
symmetric state, the two mode squeezed thermal state (TMST) , the GREE will be a minimization of a function 
on one parameter, the symplectic eigenvalue of TMST, and the state achieves the GREE is a TMST state. I 
and my coworker had calculated ^jthe minimization of relative entropy of TMST with respect to TMST as the 
upper bound of relative entropy of entanglement, now it is proved that it is just the GREE of TMST. Moreover, 
the upper bound for entanglement of formation proposed in the same paper turns out to be the entanglement 
of formation itself [TH| . So the comparison of the upper bound of EoF and RE of TMST in our former paper 
[T5] is in fact the comparison of EoF and GREE. And we had also provided coherent information and other 
entanglement measure such as logarithmic negativity in that comparison. 

I have given the method to calculate GREE for general state and the detail calculation of GREE for lxl 
gaussian state. It is expected that the method developed in this paper will be applicable to the multi-mode 
bipartite gaussian states and multipartite gaussian states. The definition of GREE need not limit to gaussian 
state. For a non gaussian continuous variable state, GREE can also be defined as the minimization of relative 
entropy of the state with respect to all separable gaussian state. But there is a deficiency that the relative 
entropy will never be zero. Never the less, the calculation involves only the first and second moments of the 
state, the necessary condition of separability on the CM of the non gaussian state was also addressed [7] , and 
the logarithmic of separable gaussian state can be treated with EM. 

Acknowledgment: Funding by the National Natural Science Foundation of China (under Grant No. 10347119) 
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